Nuclear radius deduced from proton diffraction by a black nucleus 
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We find a new method to deduce nuclear radii from proton-nucleus elastic scattering data. In this 
method a nucleus is viewed as a "black" sphere. A diffraction pattern of protons by this sphere is 
equivalent to that of the Fraunhofer diffraction by a circular hole of the same radius embedded in 
a screen. We determine the black sphere radius in such a way as to reproduce the empirical value 
of the angle of the observed first diffraction peak. It is useful to identify this radius multiplied by 
-y/3/5 with the root-mean-square matter radius of the target nucleus. For most of stable isotopes 
of masses heavier than 50, it agrees, within the error bars, with the values that were deduced in 
previous elaborate analyses from the data obtained at proton incident energies higher than ~ 800 
MeV. 
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Size of atomic nuclei, one of the most fundamental 
nuclear properties, remains to be determined precisely. 
Most popularly, the size is deduced from electron and 
proton elastic scattering off nuclei Q, 0, 0, 0, @ The 
charge radii are well determined due to our full under- 
standing of the underlying electromagnetic interactions 
HHU, while deduction of the matter radii from the 
proton-nucleus scattering data depends on the scatter- 
ing theory, which is more or less approximate in the 
sense that the nucleon-nucleon interactions involved are 
not fully understood. During the past three decades 
there have been many efforts of deducing the matter 
density distributions, which are based on various scat- 
tering theories incorporating empirical nucleon-nucleon 
scattering amplitudes, such as Glauber theory 0,13 and 
nonrelativistic and relativistic optical potential methods 
pjlsL HollIH . A systematic analysis of the data for a large 
number of nuclides, however, is still missing. In this pa- 
per we propose a method to deduce the root-mean-square 
(rms) matter radii, which is powerful enough to allow us 
to perform such a systematic analysis. This method, in 
which we assume that the target nucleus is completely 
absorptive to the incident proton and hence acts like a 
"black" sphere, is far simpler than the conventional meth- 
ods. This approximation was originally used by Placzek 
and Bethe [l^l in describing the elastic scattering of fast 
neutrons. 

The present method is useful for heavy stable nuclei 
for which the proton elastic scattering data are present, 
as we shall see. In the conventional framework to deduce 
the rms radius, one tries to reproduce empirical data for 
the differential cross section for scattering angles cover- 
ing several diffraction maxima 0,0,E1> whereas, in the 
present method, one has only to analyze the data around 
a maximum in the small angle regime. Remarkably, these 
two methods turn out to be similar in the deducibility of 
the radius. 



Elastic scattering data for more neutron-rich unsta- 
ble nuclei are expected to be provided by radioactive ion 
beam facilities, such as GSI and Radioactive Ion Beam 
Factory in RIKEN. In a possible scheme, a beam of un- 
stable nuclei, such as Ni and Sn isotopes, created in 
heavy- ion collisions is incident on proton targets, and 
the protons scattered therefrom are detected, leading to 
the measurement of differential elastic cross section. We 
expect the present method to be effective at deducing the 
rms radius of unstable nuclei from such measurement. 

We begin by regarding a target nucleus for proton elas- 
tic scattering as a black sphere of radius a. This picture 
holds when the target nucleus is completely absorptive 
to the incident proton. For high incident kinetic energy 
T p above ~ 800 MeV, the optical potential for this reac- 
tion is strongly absorptive. It can be essentially viewed 
as a superposition of the nucleon-nucleon scattering am- 
plitude. Since the imaginary part of the amplitude is 
dominant over the real part in this energy range |l4l Il5j , 
the black sphere picture is applicable to a first approxi- 
mation. We note that the black sphere picture is fairly 
successful in describing the elastic scattering of low en- 
ergy a particles 0, 0. It was also used for analyses 
of the scattering of intermediate-energy pions and low- 
energy antiprotons 0. 

Since one can regard the proton beam as a plane wave 
of momentum pLab = y (T p + m p ) 2 — m| with the pro- 
ton mass, m p , the black sphere picture can be described 
in terms of wave optics. This picture reduces to a diffrac- 
tion of the wave by a circular black disk of radius a if the 
corresponding wave optics is close to the limit of geo- 
metrical optics, i.e., a/A La b > 1, where A Lab = 2n/p Lab 
is the wave length. We will later consider the ranges of 
T p > 800 MeV and A > 50, for which a/A La b > 1 is 
satisfied. According to Babinet's principle, this diffrac- 
tion is in turn equivalent to the Fraunhofer diffraction 
by a hole of the same shape as the disk embedded in a 
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screen [17]]. The scattering amplitude for this diffraction 
in the center-of-mass (cm.) frame of the proton and the 
nucleus reads 



/(q) = ipaJ!(qa)/q, 



(1) 



where q is the momentum transfer, p is the proton mo- 
mentum in the cm. frame, and J n is the n-th order Bessel 
function. We then obtain the differential cross section 
as da/dn = |/(q)| 2 - We remark that Eq. 0) can be 
obtained from the absorptive limit of the Glauber the- 
ory in which the phase shift function is approximated by 
exp[i%(&)] = 9(b — a), where b is the impact parameter 
perpendicular to p, and /(q) is given by [Isj 

/>oo 

/(q) = % V \ bdb J Q (qa){l - exp[i X (&)]}. (2) 
Jo 

We assume that the density distribution of the black 
sphere is uniform. Then it is natural to introduce an rms 
black sphere radius, res, as 
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FIG. 1: Differential cross section calculated from the Fraun- 
hofer diffraction formula, Eq. ifT). for p- 58 Ni elastic scatter- 
ing (T p = 1047 MeV). The experimental data (crosses) are 
taken from Ref. [Til . The arrows represent the first diffrac- 
tion maximum, at which we fit the calculated peak angle to 
the empirical value. 
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In this stage, we determine a in such a way that the cm. 
scattering angle [6* c m . = 2sin (g/2p)] of the first maxi- 
mum for the Fraunhofer diffraction agrees with that mea- 
sured by proton- nucleus elastic scattering, 9m- (Here we 
define the zeroth peak as that whose angle corresponds to 
#c.m. = 0.) We remark that the diffraction patterns for 
#c.m. 3> &M are distorted by multiple scattering effects 
Pi, which are beyond the scope of the black sphere pic- 
ture. The radius, a, and the angle, 9m, are then related 
by 



2pasin(0 M /2) = 5.1356- •-. 
Combining this with Eq. 0, we may thus write 
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It is this formula that we use in the followings. As we 
shall see for heavy stable nuclei, the values of tbs that 
can be determined from Eq. JSJ agree well with the values 
of the rms matter radius, r m , deduced from elaborate 
scattering theories in previous works. For the estimate 
of r BS from the data for T p > 800 MeV and A > 50, one 
can use the approximate expression, 



r B s - 4.50 



1 GeV\ /10 deg 

PLab J \ 9m 



fm. 



(6) 



For 5S Ni and T p = 1047 MeV, the diffraction pattern 
is calculated from Eq. by setting the first peak angle 
at 6m- The result is shown in Fig. 1 together with the 
experimental data The heights of the diffraction 

maxima and minima thus calculated deviate from the 
empirical values, because the black sphere picture does 
not allow for the surface diffuseness. That is why we 



pay attention to the scattering angles of the diffraction 
maxima and minima. The relation of these scattering 
angles to nuclear sizes has been discussed by Amado et 
al. |20j in the Glauber theory. The diffraction maxima 
are more advantageous to our study than the diffraction 
minima. This is partly because statistical uncertainties 
in the yield count are much smaller near the maxima than 
near the minima and partly because theoretical predic- 
tion of the minima is rather sensitive to spin-orbit and 
Coulomb interactions between nucleons [lj. We remark 
that the behavior around the first diffraction peak can be 
well reproduced by incorporating nucleon distributions 
similar to the rectangular distribution assumed here into 
the Glauber theory (see Fig. 1 in Ref. pljl. 

It is important to note experimental uncertainty in the 
scattering angle and in the proton incident energy since 
this gives rise to the uncertainty in the estimate of res, 
together with systematic errors that are dependent on 
the way of determining the peak position. The uncer- 
tainty in the measured angle, which is due mainly to the 
absolute angle calibration, is typically of order or smaller 
than ±0.03 deg for existing data for proton elastic 
scattering off stable nuclei, while the uncertainty in the 
measured proton incident energy is typically a few MeV 

In this work we focus on the proton elastic scattering 
data for A > 50, T p > 800 MeV, and 9 c . m . > 5 deg. We 
first display the Ni results for tbs> which are obtained 
from the experimental values of 9m using Eq. Q. We 



have used the_data for T p 
and T p = 




Tn 



796 MeV 

1040 MeV |26J, and T p = 1047 MeV 0f. In collecting 
the data, we have made access to Experimental Nuclear 
Reaction Data File (EXFOR [CSISRS]) 27]. In Fig. 2, 
the values of res derived from these data are plotted 
together with the deduced values of r m . When the values 
of r m are not explicitly given in the literatures, we obtain 
them from = (Z/A)r^ + [(A - Z)/A]r*, where Z is 
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FIG. 2: (Color) The values of rss (squares) for Ni isotopes. 
For the error bars, see text. The values of T p (in MeV) are 
specified above these error bars. For comparison we also plot 
the results for r m derived in the following references: For 
A = 58j_in Refs. El |H |gl IH |H (T p = 796 MeV) and in 
Refs. E l ll ljpl (T p = 1047 MeV); for A = 60 and 62, in 
Refs. (HlSll; for A = 64, in Refs. (T p = 796 MeV) 

and in Refs. TUlH |H (T p = 1047 MeV). The crosses (x) 
denote the rms matter radii of the point nucleon distributions, 
and the circles (o) denote those folded with the nucleon form 
factor. Only r B s is plotted for A = 58 and T p = 1040 MeV, 
because no r m is available for this case. 



the charge number, and r p and r n are the rms radii of 
the proton and neutron distributions. The values of 9m 
are determined from the scattering angle that gives the 
maximum value of the cross section among discrete data 
near the first diffraction maximum. The error bars of 
res in the plot, which are of order or even greater than 
±0.05 fm, are determined from the half width between 
the neighboring angles measured. Uncertainties of r B s 
associated with the absolute angle calibration and the 
measurement of T p , which are smaller than ±0.03 fm and 
±0.02 fm, respectively, are not taken into account. When 
the error of r m is not available, we evaluate it from the 
errors of r n and r p given in the literatures. 

We find from Fig. 2 that for all the isotopes, the val- 
ues of res agree with those of r m within the error bars 
except for an only case. This result is remarkable since 
there is no fitting parameter other than a in the black 
sphere picture in contrast to the past elaborate analyses. 
We note that it is premature to ask whether we have to 
regard r B s as an rms radius of the point nucleon distribu- 
tion or as an rms matter radius folded with the nucleon 
form factor. In fact, the present estimate of r B s contains 
errors larger than the difference between these two radii, 
which is typically ~ 0.07 fm. 

We next show the relations between r B s and r m in 
Fig. 3, which are constructed from systematic data for 
various stable nuclides of mass number larger than 50. 
In deriving r B s we have used the data for 54 Fe (T p = 796 
MeV) |jt 90 Zr (T p = 800 MeV) |H, 90 Zr (T p = 1000 
MeV) 36], 116 ' 124 Sn (JL = 800 MeV) H2, and 208 Pb 
(T p = 800 MeV) (HQ, 208 Pb (T p = 1000 MeV) [36|, 
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FIG. O'. TBS VS. Vm for stable nuclei of masses above 50. For Ni 
isotopes the values of r m are taken from the same references 
as cited in Fig. 2; for 54 Fe, from Ref. El; for 90 Zr, from Refs. 
PlSlilrfor 116 > 124 Sn^ from Refs. |M |H H2; for 208 Pb, 
from Refs. [II |H US IH IH |H . The definition of the 
crosses (x) and the circles (o) is the same as in Fig. 2. The 
dotted line represents res = r m . Inset: tbs vs. r m for Ni 
isotopes. 



We do not include the data of 208 Pb (T p = 1040 MeV) 
|26j | in this analysis, because the first peak position is not 
clear. The data for Ni isotopes are the same as used in 
Fig. 2. We can see that the values of res and r m including 
the error bars are mostly on the line of res = f m . We 
thus find that res provides a good measure of the rms 
matter radius. If elastic scattering data are obtained in 
a much finer manner, the main uncertainty in res would 
arise from the absolute angle calibration and possibly 
the measurement of T p . In any case one could nicely 
determine the isotope dependence of res if the relative 
peak angles between isotopes are accurately measured for 
the same proton beam. 

For nuclides for which elastic scattering data are avail- 
able but no r m , we list the values of tbs for reference. For 
90 Zr (T p = 800 MeV) H3, we obtain r BS = 4.22 ± 0.10 
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for 92 Zr (T p 
120 Sn (T p = 
144 Sm (Tp -- 
154 Sm (T p -- 
176 Yb (T p = 
208 Pb (T p - 



90 Zr (Tp = 800 MeV 
800 MeV) la 




|3J, r BS = 4.21 ±0.11 fm; 
, r BS = 4.21 ± 0.07 fm; for 
r BS = 4.62 ± 0.10 fm; for 
r BS = 4.92 ± 0.27 fm; for 
r BS = 5.24 ± 0.09 fm; for 
r BS = 5.47 ± 0.10 fm; for 
r BS = 5.54±0.29 fm. Among 
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these nuclei, 154 Sm and 176 Yb are deformed in the ground 
state. Since cross sections are generally measured for pro- 
ton elastic scattering from randomly oriented nuclei, the 
disk radius, a, obtained from the data lies between the 
lengths of the semimajor and semiminor axes of the de- 
formed nuclei. As long as the degree of deformation is 
small, therefore, r B s is expected to give a good measure 
of the rms matter radii of the deformed nuclei. 

In summary we have performed a systematic analy- 
sis, based on the Fraunhofer diffraction off a black disk, 
of the existing data for proton elastic scattering off sta- 
ble nuclei with masses larger than 50 at proton incident 
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energies above ~ 800 MeV. This analysis allows us to 
make systematic estimates of the rms matter radii. The 
present method works even for the data in a range of q 
covering only the first diffraction maximum. Such cases 
would occur for neutron-rich unstable nuclei. By com- 
bining the present result with a comprehensive table for 
the rms charge radii 0, 0] , we can estimate neutron skin 
thickness for various nuclides on equal footing. This is 
an important step towards understanding of the isospin- 
dependent bulk and surface properties of nuclear matter 

In the near future, an experiment that provides dif- 
ferential cross sections of proton elastic scattering off Ni 
unstable isotopes will be performed in GSI. In this exper- 
iment the projectile will be a radioactive ion beam having 
energy per nucleon of about 400 MeV, and the scattering 
angles to be measured will contain the first diffraction 
maximum. In applying the present prescription to esti- 
mate the rms matter radius to such measurements, its 
validity, which has been confirmed here for proton inci- 
dent energies above ~ 800 MeV, needs to be examined 
at relatively low bombarding energies. Investigation of 
how res depends on the bombarding energy is now in 



progress |4l|. Once the values of res are accumulated 
for various nuclides including neutron-rich heavy nuclei, 
it is expected that the density dependence of the symme- 
try energy near normal nuclear density will be clarified 
|42|. This expectation is also suggested by the work [2l| 
that pointed out the relation of 9m with the density de- 
pendence of the symmetry energy. We note, however, 
that the estimate of res depends strongly on how sharp 
the energy distribution of a radioactive ion beam will be. 

In order to deduce the surface diffuseness in addition 
to the rms radius, one has to reproduce the overall be- 
havior of the elastic scattering differential cross sections, 
as shown in Refs. 0,0,0]. The tails of the distributions, 
which are hard to deduce from the differential cross sec- 
tions, might affect the prediction of reaction cross sec- 
tions 43]. If reaction cross sections are measured for 
various heavy nuclides, the "halo" structure of unstable 
nuclei and the tails for stable nuclei will become clearer. 
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